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Abstract 

Precise calculations of the coexistence-curve diameters of a hard-core square-well (HCSW) fluid and the restricted primitive 
model (RPM) electrolyte exhibit marked deviations from rectilinear behavior. The HCSW diameter displays a |t | singularity 
that sets in sharply for \t\ = \T-Tc\/Tc < 10"^; this compares favorably with extensive data for SFg, also reflected in C2H6, 
N2, etc. By contrast, the curvature of the RPM diameter varies slowly over a wide range \t\ < 0.1; this behavior mirrors 
observations for liquid alkali metals, specifically Rb and Cs. Amplitudes for the leading singular terms can be estimated 
numerically but their values cannot be taken literally. 



1. Introduction 

Coexistence curves of gas-liquid phase separation 
in single-component fluids have been investigated ex- 
tensively for more than a century via both experiment 
and theory. In particular, the van der Waals (vdW) 
equation of state has been successful in describing 
the experimental data at a semiquantitative mean- 
field level [1]; however, the shape of the coexistence 
curve predicted by the vdW equation of state dif- 
fers significantly from observation. Nevertheless, one 
aspect of the vdW equation of state is that the diam- 
eter of the coexistence curve, Pdmm{T), namely the 
mid-points of the coexistence-curve boundary in the 
density-temperature plane, exhibits a linear behavior 
with temperature, T, as the critical point, Tc, is ap- 
proached; this is the Law of the Rectilinear Diameter, 
and is typically observed experimentally to be valid 
with a fairly high degree of precision. 



Several models for fluids, however, have been ad- 
vanced that violate the law [2,3,4]. These models sug- 
gest that the diameters of real fluids should exhibit an 
entropy-like singularity varying as near criti- 

cality where t = {T - Tc)/Tc and a ~ 0.109 is the 
critical exponent for the specific heat; the slope of the 
diameter was thus expected to diverge at criticality in 
the same fashion as the specific heat. The traditionally 
accepted scaling theory of the critical region, which 
incorporates mixing between the temperature T and 
the chemical potential fi in the scaling fields (but does 
not include the pressure), also generates this entropy- 
like singularity [5,6]. However, the recent discovery 
of a Yang- Yang anomaly [7,8] in fluid criticaUty turns 
out to entail a further singular term in the diameter 
proportional to where (3 ~ 0.326 is the criti- 

cal exponent for the order parameter [7,9]. And such 
a term also arises in other classes of exactly solu- 
ble model [7,10]. In fluid criticality, one normally has 
2/3 ~ 0.65 < 1— a so that the \t\'^^ term dominates 
the contribution when t ^ 0— . 



Corresponding Author: xpectnil@ipst.umd.edu 



Preprint submitted to Elsevier Science 



2 Februmy 2008 



Deviations from the Law of the Rectilinear Diame- 
ter have indeed been observed in several careful exper- 
iments on various fluids [11,12,13,14,15,16]. Among 
these, the seminal study of SFe by Weiner et al. [12] 
more than 30 years ago reported a surprisingly sharp 
and abrupt singular downturn of the diameter very 
close to Tc- More specifically, their data reveal a fairly 
linear behavior (as predicted by the vdW equation) 
over a wide range of temperature up to t ~ —10^'^; 
but then the diameter drops sharply to the observed 
critical density, pc- see Figs. 1(a) and 2(a), below. This 
sudden drop seems to reflect the anticipated 
singularity, provided one may assume that the Yang- 
Yang anomaly is negligible in this system. 

Subsequent experiments on other single-component 
fluids, specifically, C2H6 and N2 have revealed similar 
sharp singularities, while, although less definitively, 
they have also been observed in C2H4 and Ne [15,16]. 
To our knowledge, however, no theoretical or numeri- 
cal calculations for model fluids have yielded such be- 
havior. One might thus reasonably harbor some doubts 
about the SFe and similar data: see also a potential ex- 
perimental concern [17]. Nevertheless, Ley-Koo and 
Green [18] analyzed these data for the diameter us- 
ing the traditional restricted scaling ansatz with no 
pressure mixing (i.e., ji = = 0: see below). They 
found that the singular \ t\^~°' term with corresponding 
higher order corrections sufficed to describe the data 
for |f| < 10~^ reasonably well with a ~ 0.15 (which 
is rather larger than the currently accepted value a ~ 
0.109). 

On the other hand, the experimentally observed co- 
existence curves for liquid alkali metals, such as K, 
Rb, Cs, etc., display quite different behavior from that 
of SFg [13,14]. These metals exhibit marked devia- 
tions from the Law of the Rectilinear Diameter with 
a smooth downward curvature extending over quite a 
wide range of temperature even below T ^ 0.9 T^, but 
they do not show a sharply dropping singularity near 
criticality as seen for SFg: see Fig. 3, below. Clearly, 
in Ught of the varied experimental data it is of inter- 
est to see whether specific model fluids, that perhaps 
capture some aspects of real systems, might produce 
such disparate types of behavior. 

Until recently, calculations for continuum (or off- 
lattice) model fluids, mostly computer simulations, 
have failed to observe singular behavior in the diame- 
ter such as seen in the experiment on SFg or even the 



extended curvature of the alkaU metals. However, a 
newly developed scaling algorithm has enabled us to 
obtain precise coexistence curves for continuum mod- 
els on the basis of grand canonical Monte Carlo (MC) 
simulations [19,20]. Equipped with suitably designed 
finite- size scaling techniques [19,21], the unbiased al- 
gorithm has generated [22] the coexistence-curve di- 
ameters up to 1 part in 10'^ or 10'* of Tc for a hard-core 
square-well (HCSW) fluid with an attractive-well of 
range 1.5(T (where a is the hard-core diameter of the 
particles) and for the restricted primitive model (RPM) 
electrolyte — equisized hard spheres half with charge 
+q and half with —q. As we report here, the simulation 
results for the diameter for the HCSW fluid turn out 
to resemble surprisingly closely the experimental data 
for SFg, exhibiting a similar sharp drop very close to 
Tc! By contrast, the estimated diameter for the RPM 
is found to mirror the observed behavior of the liq- 
uid metals: specifically, /9diam(2^) curves slowly down- 
wards towards pc as Tc is approached with no visibly 
obvious singularity or rapidly growing slope. In both 
the RPM and the metals the role of the Coulombic in- 
teractions is felt to be paramount; so this difference 
in behavior, which seems also to be correlated with a 
large value of the Yang- Yang ratio, "Jifj,, was, perhaps, 
to be expected. 

Our goal here is to compare quantitatively the cal- 
culated diameters of flie HCSW fluid and the RPM to 
the experimental data for SFg and the liquid metals 
[12,13,14]. The comparison is instructive in that one 
can gauge how well these simple models capture some 
essential features of real fluids and, at the same time, 
throw hght on various experimental [17] and theoret- 
ical [23] proposals. 



2. Theoretical background 

For completeness and perspective we sketch here 
some of the underlying theory. The thermodynamic 
potential of a one-component fluid near the critical 
point {Pc,Tc,i^c) can be written in general scaUng 
form with the leading correction as [7,9] 

p « Q\i\^-''W±iUh/\i\^, Ui\if), (1) 

with ± for t ^ 0, while, to Unear order, the relevant 
scaling fields are 
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p = p - kot - lojj. + ■ 
i = t-hjl- jip + ■ 
h = - kit - j2P + ■ 



(2) 
(3) 
(4) 



in which we have introduced the dimensionless critical 
deviations 



keTc 



P-Pc 
Pck^Tc 



(5) 



while in (1) the gap exponent is A = 2 — a — (3 = 
(3 + 7 where, in standard notation, 7 is the exponentfor 
the divergence of the compressibihty/susceptibility, 
while 6 is the leading correction-to-scahng exponent. 
Together with the exponents a, fi, ■ ■ •, etc., the scal- 
ing functions W±{x,X4) are universal. On the other 
hand, the amplitudes Q, U, and U4 are nonuniversal, 
i.e., system-dependent as are the mixing coefficients 
ji,j2,---,h- In this expression we have neglected 
higher order corrections and nonUnear mixing terms 
in the scaling fields [9]. 

The scahng ansatz (l)-(4) accounts for the Yang- 
Yang anomaly: explicitly, when j2 is nonvanishing, the 
second temperature derivative of the chemical poten- 
tial on the phase boundary, ij.'^{T) = {dP^^^/dT'^), di- 
verges at criticaUty like the specific heat, i.e., as Cy ^ 
The Yang- Yang ratio is then given by [7,9] 

3?^ = - ^lim_ T,,';iT)/Cv{T) = 32/ {j2 - 1). (6) 

The two sides of the coexistence curve, p = 
p'^{T) > Pc and p^{T) < pc, follow directly from 
the scaling ansatz (1). The results are [9] 



Ap{T)^[p+{T)-p-{T)]/p, 

= B\tf{l + be\tf + ■■■), 



(7) 



A/5d(r) : 



p+{T)+p-{T) 
2pc 



= A20\t\'"^ + Ai-a,\t\^-" +Ait 

+A'^^\tff'+' + Ae\t-''+' + ---, 



(8) 



where the amplitude B, and coefficients bg, A213, 
etc., can be expressed expUcitly in terms of universal 
derivatives of W± (x, 2:4) and the nonuniversal values 
of Q, [/, Ui and the mixing coefficients in (2)-(4). 
Note that is proportional to the pressure-mixing 
coefficient 32 [9], while Ai-a is proportional to 



'1 + ii- For the {d = 3) Ising universaUty class, which 
describes standard fluid criticality, one has a ~ 0.109 
and f} ~ 0.326, as mentioned, while 9 ~ 0.52 [24]. 



3. Diameters of SFe and the HCSW fluid 

The coexistence-curve diameter of SFg measured 
by Weiner et al. [12] is presented in Fig. 1 together, 
in order to make a comparison, with the simulation 
results for the HCSW fluid [25,22]. Both observation 
and calculation span almost three orders of magnitude 
in temperature reaching up to within a few parts in 
10"' of Tc. It is striking, first, to see that the HCSW 
fluid, the most basic continuum model that undergoes 
gas-liquid phase separation, exhibits a behavior that 
so closely matches that for SFe. Indeed, the diameters 
of both systems show essentially linear behavior for 
\t\ > 0.005, as may be gauged by the pair of parallel 
straight, dot-dashed lines that embrace the data in both 
Figs. 1 (a) and (b). Note, however, that extrapolation 
of these rectilinear regimes to criticality yields density 
offsets, Apc, above the true value of pc amounting to 
about 0.53% for SFg and 0.23% for flie HCSW fluid. 
Then, the diameters for both systems display unex- 
pectedly sharp downwards departures from Unearity 
that set in around t = —(3-5) x 10~^ and suggest a 
divergent derivative at Tc. As far as we are aware this 
is the first time that such seemingly singular behavior 
in the diameter of a model fluid has been identified 
via computer simulations. 

In order to reveal the putative singularities in more 
detail. Fig. 2 presents magnified plots of the diame- 
ters over the much smaUer temperature interval \t\ < 
0.005. Although the diameter of flie HCSW fluid ex- 
hibits a somewhat sharper drop, the quantitative simi- 
larity in this wowuniversal aspect of the critical behav- 
ior is surprising granted the crude nature of the HCSW 
fluid regarded as a model of a real system undergoing 
gas-liquid phase separation. The close correspondence 
between the observations and the theory suggests that 
the issue raised by Moldover and Gammon [17] (and 
cited in [23]; but see also [16]) regarding the possi- 
ble influence of wetting layers on the observations of 
Weiner et al. should not be regarded as a serious con- 
cern. 

In order to understand the contributions of the sin- 
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Fig. 1. Plots of the reduced coexistence-curve diameters, 
Apj,(T) s (Pdiam/Pc) - 1 vs t = (T - Tc)/Tc for (a) SFg, as 
observed by Weiner et al. [12], and (b) the hard-core square-well 
(HCSW) fluid with interaction range 1.5a, as determined via 
Monte Carlo simulations [22]. For clarity, none of the data points 
closer to Tc than 2-3 part in lO'' are shown in either plot while 
most of the calculated data points for \t\ < 0.01 have been omit- 
ted in (b): these data can he examined in Fig. 2. The solid curves 
represent fits to all the data as detailed under (ii) in Table 1. 
Conversely, the dot-dashed lines that narrowly encompass the data 
points serve to demonstrate that the Law of the RectiUnear Diam- 
eter is well obeyed for t > 0.005. Note, however, that this linear 
behavior extrapolates to offsets. Ape, above pc , as marked by 
the open arrowheads on the ordinate scales. 

gular terms to the diameter, we have examined fits of 
the data to (8) subject to various conditions. First, we 
may neglect the term, since its amplitude A'^/^ 

is proportional to j2 and will probably be of small if 
any significance in light of the estimated and likely 
value of the Yang- Yang ratio, 51^ [7,22]. Next, on 
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Fig. 2. The reduced coexistence-curve diameters, Apjj(r), close 
to Tc for (a) SFe [12] and (b) for the HCSW fluid [22]. The 
SFe data presented here [and in Fig. 1(a)] are reproduced from 
Figs. 1 and 2 of Weiner et al. [12] (with neglect of the two data 
points closest to Tc that were subject to systematic displacements 
entailed in linear extrapolations for the gravitational corrections 
implemented). The various curves represent fits to the data based 
on the expression (8) with the coefficients given in Table 1. Note 
that the dot-dashed plot (iii) in (a) is almost identical to the soUd 
curve, the differences being visible only for t > 0.003. 

ignoring, initially, the leading term varying as \t\'^^, 
that originates from the pressure mixing (and so is 
also proportional to 3?^), we find least-squares fits that 
seem to be significantly less than satisfactory: see the 
dashed fines labeled (i) in Fig. 2. The associated val- 
ues of the coefficients for SFg are about twice those 
for the HCSW fluid: see lines (i) in Table 1 . Of course, 
this might well have been anticipated from Fig. 2 in 
which the vertical scale for SFg is twice that for the 
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Table 1 

Parameters for fitting tiie diameters of SFg [12] and the HCSW 
fluid [22] to Eqns. (8) and (9). 







A213 




Ai 


Ae 




(i) 


SFs 




6.365 


-10.13 


8.080 






HCSW 




2.897 


-4.532 


3.884 




(ii) 


SFe 


1.124 


-9.042 


11.37 


-3.354 






HCSW 


0.965 


-9.818 


12.97 


-5.092 




(iii) 


SFe 


1.086 


-7.990 


9.770 




3.318 




HCSW 


1.143 


-11.15 


14.38 




12.02 



HCSW fluid. Nevertheless, the fact that the ratio of 
corresponding coefficients lies within the rather nar- 
row range 2. 16 ± 0.08 provides an interesting measure 
of the similarity of the two data sets. 

On the other hand, if one also allows the leading, 
pressure mixing term in (8), the fits are markedly im- 
proved: see the solid curves in Figs. 1 and 2. Further- 
more, the values of the three leading fitted coefficients 
are now surprisingly close for the two systems; in- 
deed, only the last amplitude, Ae (which pertains to 
the leading corrections to scaUng) differs by more than 
17%: see (ii) in Table 1. 

To gain some insight as to what might be reasonable 
confidence limits to attach to these results, we have 
fitted the data to the generahzed Pade approximant 
form 

1 + ae\t\'' 

where the \t\^ term in the denominator accounts for 
infinitely many higher order contributions (including 
the neglected term). The results of this fit are 

presented in lines (iii) in Table 1 and as dot-dashed 
curves in Fig. 2. As seen in the figure, the differences 
between the fits (ii) and (iii) are negligible. Further- 
more, the values of the three leading coefficients have 
changed by no more than 20%. 

Now, according to the scaling ansatz (l)-(4), the 
leading amplitude in (8) can be expressed as [9] 

A20 = ^i,B\ (10) 

where IR^ is the Yang- Yang ratio defined as in (6). This 
ratio and the amplitude B have been estimated rather 
precisely for the HCSW fluid [19,22]. Using the values 
~ -0.04 and B ~ 2.0 then yields A2/3 ~ -0.16. 



This value, however, is not only an order of magni- 
tude smaller than the fitted values in Table 1 , but also 
has the opposite sign! This indicates unequivocally 
that the fitted coefficients are unlikely to represent the 
true amplitudes for the diameter Possible reasons are, 
first, that since the contributing terms in (8) have ex- 
ponents that are closely spaced numerically, they in- 
terfere strongly in the fitting procedure, resulting in 
effective rather than realistic amplitudes. Second, one 
knows that there are higher order terms which have 
been totally ignored in the fitting but that might play 
a crucial role. One such is the leading odd correction- 
to-scaling contribution varying as [9]; others 
will arise from the nonlinear contributions to the scal- 
ing fields (2)-(4) [9]. Finally, if one fixes at the 
expected value —0.16 and considers the approximant 
(9) with an extra term in the numerator, 
one finds that the diameter of the HCSW fluid can be 
described just as well with the coefficients Ai^a — 
206, Ai ~ -377, A'g ~ 1144 and ae ~ 1872 as with 
the sets (ii) and (iii) in Table 1 . From this lesson, we 
learn that beyond demonstrating the close similarity 
of the two systems, one must not attach much signifi- 
cance to the actual coefficient values for SFg and the 
HCSW fluid in Table 1 . In order to obtain reliable es- 
timates for even the two leading singular amplitudes, 
A2P and Ai-a > one would require measurements of 
the diameter with markedly higher precision and still 
closer to Tc. 

Finally, we mention that Goldstein and coworkers 
[15,16] examined the trends with molecular proper- 
ties of both the slopes of the coexistence-curve diam- 
eters and the strengths of the singularity for 'normal 
fluids,' [26] such as SFg, C2H6, N2, etc. They argued 
that a crucial role is played by triplet or three-body 
intermolecular interactions that, in turn, may be quan- 
tified by the molecular polarizabilities [15,16]. (Ob- 
served diameter slopes had previously been correlated 
with quantum-mechanical effects as embodied in the 
De Boer parameter A* = h/ Vmu^e, where m, a and 
e represent the molecular mass, collision diameter and 
attractive well depth, respectively [27].) In contradic- 
tion, however, later measurements for xenon [28], for 
which the polarizability product, UpPc [15,16] is un- 
usually high, reveal a diameter of large slope but un- 
detectably small deviations from rectilinearity close 
to Tc. Furthermore, Singh and Pitzer [26] studied a 
larger array of normal fluids and concluded generaUy 
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Table 2 

Parameters for fitting the liquid and vapor densities of the alkali 
metals [13,29] and the RPM [19,22] to Eqns. (7) and (8). 





Mil 




Ai 


Ae 


B 


be 


Rb 


-5.041 


54.65 


62.73 


17.21 


1.908 


0.556 


Cs 


-0.232 


2.206 


-0.611 


-0.911 


1.927 


0.436 


RPM 


-3.120 


45.85 


47.09 


9.564 


6.932 


0.312 



that the slope of the diameter does not correlate with 
three-body forces, but, rather, is primarily determined 
by the shape of the two-body potential. That, in turn, 
relates closely to the acentric factor defined via the 
slope of the vapor pressure curve at T = 0.7Tc [26]. 
In as far as there are no many-body interactions in 
HCSW models, our results for the singularity are in 
accord with this conclusion. 



4. Diameters of liquid metals and the RPM 

The coexistence curves of the liquid alkali metals 
[13,14] are found to exhibit a large degree of asym- 
metry and so are quite distinct in their overall shapes 
from those of SFg and other normal fluids and from the 
HCSW model. See Fig. 3 for Cs and Rb and Hensel 
et al. [14] for Na and K. This behavior can be com- 
pared to that of file RPM [19,22,30,31] which is even 
more strongly asymmetric. The relative similarity of 
the coexistence curves of the Uquid metals and the 
RPM would appear to be instructive since the long- 
range Coulombic interactions play a dominant role in 
both types of system. 

In order to obtain the diameters of the liquid met- 
als, we have performed a global least-squares fit to 
the coexisting densities, p'^iT), on the basis of the 
expansions (7) and (8). The values of Tc and pc used 
in our fits for Cs and Rb are those stated by Hensel 
and coworkers [13] (But see further below). The re- 
sulting amplitudes and coefficients are presented in 
Table 2 together with the corresponding RPM values. 
The coefficients • • •, Ag for the diameter of Rb 
prove quite comparable to those of the RPM; however, 
the leading amplitude, B, for the density discontinu- 
ity, Ap{T), of the RPM is larger by a factor close to 
3.6 which simply reflects the markedly greater asym- 
metry of the model relative to the liquid metals and, 
concomitantly, also entails a more steeply sloping di- 
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Fig. 3. The coexisting densities of Cs and Rb as reported by Hensel 
and coworkers [13] who invoked further data for Cs (crosses) 
from Stone et al. [29]. The solid curves portray fits to (7) and (8) 
while the dashed curves represent the corresponding diameters: 
see Table 2 for the coefficients employed. Note that the tabulated 
coexistence curves of these liquid alkaU metals [13] were obtained 
by extending isochores down to the previously observed vapor 
pressure curves. Since, as a result, the liquid and vapor densities 
are available only at differing temperatures, the diameters cannot 
be determined directly and fitted only to (8). 



ameter [22,30]. 

The diameters of the liquid metals determined by 
these procedures and of the RPM (as found in [22]) are 
presented in Fig. 4. It is striking that, apart from the 
anticipated multiplication factor close to 3, the diam- 
eter of the RPM is almost identical to that of Rb. Both 
diameters exhibit smoothly varying curvature down 
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Fig. 4. The reduced coexistence-curve diameters, Ap^{T), of Cs 
(dashed), Rb (sohd curve), and the RPM (solid circles). The plots 
for Cs and Rb, derived from [13] by fitting to (7) and (8), have 
been amplified by a factor 3 to facilitate comparison with the 
calculations for the RPM for which the uncertainties are judged 
to be no larger than twice the symbol size [22]. The dot-dashed 
curve for Cs is found by changing slightly (within acceptable 
uncertainty Umits) the stated values for Tc and pc [13]: see text. 

to t ~ —0.2 or lower and differ significantly from 
those of SFe and the HCSW fluid. The diameter of Cs 
exhibits a less pronounced curvature, but still shows 
distinctive departures from rectilinear behavior as, in- 
deed, originally stressed by Hensel and coworkers. It 
should be noted, however, that if one allows small 
deviations from the quoted critical parameters [13], 
namely, say ~ 0.02% of (or only < 0.4 K) and of 
1% in (which seems allowable in light of the tab- 
ulated data), the revised fits reproduce curvatures as 
large as those for Rb: see the dot-dashed plot in Fig. 4. 

Just as for the HCSW fluid, the leading amphtude, 
A213, of the diameter for the RPM can be estimated 
via (10) by using the value ~ 0.26 [19,22] and 
the amplitude B given in Table 2: this yields ^2/3 — 
12.5 which is four times the magnitude of the fitted 
coefficient in Table 2, and, more disturbingly, has the 
opposite sign. Once again, this demonstrates that fitted 
amplitudes, such as recorded in Table 2, may not be 
regarded as reaUstic. At best, owing to the interplay of 
closely spaced exponent values and contributions from 
higher order corrections, the data in Table 2 must be 
regarded as effective values that, collectively, describe 
the observations reliably. This is further supported by 



the fact that fixing A2/3 to the theoretical value 12.5 

and using the approximant (9) yields a satisfactory fit 
with the coefficients Ai-a — —199, Ai ~ —327 and 
ae ~ 35. 

Two further remarks are called for. First, it must be 
remembered that the sudden downward turns of the 
diameters of SFg and the HCSW fluid set in around 
t = —(2-3) X 10~^ and are clearly developed only 
for \t\ < (4-5) X 10"^: see Figs. 1 and 2. Owing to 
experimental difficulties, however, reliable data for the 
alkali metals have been obtained only for \t\ > 10~^ 
[13]. It is possible, therefore, that further observations 
approaching Tc by another decade or more would also 
uncover a sharp drop in the liquid metal diameters 
very close to Tc {in addition to the extended curved 
region further from criticality). 

Second, the experiments of Hensel and coworkers 
were associated with an interesting theoretical sug- 
gestion of Goldstein and Ashcroft [23]. They argued 
that the characteristic curvature of the coexistence- 
diameters of the liquid metals (in contrast to the 
largely rectilinear behavior observed for typical in- 
sulating or normal fluids) "arises from correspond- 
ingly strong thermodynamic-state dependence of the 
screened ion-ion interactions" in such systems. In ad- 
dition, they suggested an association with the crossing 
of the metal-insulator transition. As regards this sec- 
ond proposal, we believe the weight of evidence from 
various studies now suggests that the metal-insulator 
crossover (which is not a sharp transition at nonzero 
temperatures) has essentiaUy no correlations with 
thermodynamic properties beyond residing in the wide 
range of fluid densities accessible near and above any 
gas-hquid critical point: see, e.g., studies of the phase 
behavior of Uquid mercury through criticaUty [32]. 

However, the close similarity of the coexistence di- 
ameters of the liquid metals and the restricted primitive 
model (RPM) electrolyte casts a fresh Ught on the first 
suggestion. In as far as the RPM is a perfectly standard 
statistical mechanical model with well-defined /zxeii, 
i.e., state-independent intCTaction potentials, albeit of 
long range, the focus of Goldstein and Ashcroft on 
state-dependent potentials per se might seem inappro- 
priate or, at least, unnecessary. Certainly, neither quan- 
tum mechanics, strongly delocaUzed electronic states, 
or rapid changes from localized tight-binding to itin- 
erant band states, nor the associated electron-ion mass 
and interaction disparities can be directly relevant! 
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On the other hand, if one wishes to view a classi- 
cal ionic systems as an uncharged fluid with effective 
"screened interactions" one is naturally led to con- 
sider pair potentials of Debye-Hiickel form. The in- 
teraction range is then set by the charge screening 
length, £^z{T,p), which is certainly state-dependent, 
approaching the Debye value, oc {T/ pY/"^ at low 
density and high temperatures [33]. However, Gold- 
stein and Ashcroft [23] characterize such screening in 
"typical electrolyte solutions" as only "weakly state 
dependent" in contrast to the "vast qualitative change 
in the form of the interaction which occurs at the 
metal-insulator transition." Although theoretical argu- 
ments [34] indicate that, even in the fully charge- 
symmetric RPM, weak energy-like singularities ap- 
pear in various charge correlation lengths near criti- 
cality, it is hard, in our opinion, to view these changes, 
even if relatively rapid, as driving the curvature and 
seeming singular behavior of the coexistence-diameter 
in the RPM or, mutatis mutandis, in the liquid metals. 

Nevertheless, the Coulomb dominated behavior of 
the RPM and of the liquid metals does appear to be 
directly associated with the coexistence-curve asym- 
metries and, it even seems plausible, with relatively 
large magnitudes of the Yang- Yang ratio Dl^. Further 
experiments, simulations, and calculations to investi- 
gate these surmises would surely be of interest. 



5. Summary 

In summary, the experimentally observed singular 
behavior or, more precisely, the deviations from recti- 
linear variation of the coexistence-curve diameters of 
both SFg and liquid alkali metals have been compared 
with rather precise numerical calculations for a hard- 
core square-well (HCSW) fluid and for the restricted 
primitive model (RPM) electrolyte [22]. The close 
similarity between the accurately measured diameter 
of SFg [ 1 2] and that of the HCSW fluid is remarkable; 
the sharp singularity setting in very close to is strik- 
ing and resembles the entropy-like singularity pro- 
posed long ago even though fits to the data are unable 
to confirm the specific analytic expectations unequivo- 
cally. By contrast, the markedly curved diameters dis- 
covered for the Uquid alkali metals [13] exhibit very 
different behavior from those of SFg and the HCSW 



fluid; it transpires, instead, that they closely reflect the 
calculated behavior for the RPM. While it seems likely 
that the differences from the non-conducting systems 
are intimately associated with the long-range Coulomb 
interactions, proposals [23] emphasizing the delocal- 
ized, quantum-mechanical electronic states normally 
used to describe liquid metals, do not seem pertinent. 

Finally, let us mention that the rapidly dropping sin- 
gular behavior seen in the coexistence-curve diameters 
for the HCSW fluid and for SFe and various single- 
component fluids has not so far been observed in sim- 
ulations of other models. A similar sharp near-critical 
behavior has been reported in the diameter of a binary- 
liquid system [35] and, indeed, the same spectrum of 
singularities is to be expected in general. However, 
because of the extra thermodynamic degree of free- 
dom in a binary system, a further set of mixing co- 
efficients for the second chemical potential must be 
invoked. As a result both the theoretical scaling anal- 
ysis and the quantitative elucidation of experimental 
observations face extra challenges. Nevertheless, fur- 
ther experimental studies of such systems and, indeed, 
corresponding calculations could be instructive. 
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